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Abstract: We highlight a few recent results on the effect of the diffusion process in deterministic area preserving maps 
with noncompact phase space, namely the standard map. In more detail, we focus on the anomalous diffusion arising due 
to the accelerator modes, i.e. resonant-like features of the phase space which are transported in a super-diffusive (ballistic) 
manner. Their presence affects also the trajectories lying in the immediate neighborhood resulting in anomalous (non- 
Gaussian) diffusion. We aim to shed some light on these special properties of the phase space by utilizing the diffusion 
exponent (rate of diffusion) and the momentum distribution in terms of the Levy stable distributions, with the goal to 
reach an understanding of the global behaviour. To this end we consider a rather large interval for the kick parameter 
(0 < K < 70) of the standard map where accelerator modes of different periodicity exist. 


1. INTRODUCTION 

In this paper we present a small review of our recent 
findings on the diffusion properties in area preserving 
maps, exemplified by the standard map [1], This work 
has been linked to the extensive studies in [2,3] of the 
quantum kicked rotator (see e.g. [4]), where the classical 
diffusion process also plays a significant role in the quan¬ 
tum semiclassical level, where the (exponential) quantum 
(or dynamical) localization is observed and studied. We 
here present a few results on how the presence of accel¬ 
erator modes affects the dynamics of the classical phase 
space, and in more detail its diffusion rates and transport. 
We pursue this goal by studying the diffusion process of 
ensembles of initial conditions for a rather large value- 
interval of the control (kick) parameter (0 < K < 70). 
Within this interval the phase space of the model shows 
a very rich variety of dynamics, namely from the almost 
entirely regular phase space to gradually more and more 
chaotic dynamics until there are practically no stable re¬ 
gions anymore. In this picture there are some resonance 
features, the so-called accelerator modes, that may affect 
the diffusion of trajectories rather drastically. 

Here, we are trying to shed some light on these effects 
by calculating the diffusion exponent p, the correspond¬ 
ing diffusion constant 1) f , and the parameters of the rele¬ 
vant Levy stable distribution in the case of non-Gaussian 
(anomalous) diffusion for different K parameter values. 
Comparing to previous studies, here and mainly in [5] 
we do not restrict ourselves to local effects and isolated 
accelerator modes but we are interested in a rather more 
global point of view. Namely, on how these features may 
influence the transport in their neighborhood. That is the 
case when one, for example, considers mixed ensembles 
of initial conditions, e.g. ensembles that contain accelera¬ 
tor modes together with chaotic trajectories around them. 


Let us also stress here that the time scales we are inter¬ 
ested in in this particular study are relatively short and 
we are not interested at this point on the asymptotic dif¬ 
fusion rates. This has to do with our initial motivation in 
understanding the quantum effects associated with such 
features, namely the quantum localization properties (see 
[2,3] for more details). 

The paper is structured as follows: In Sec. 2 we give 
brief description of the model (standard map), while in 
Sec. 3 we present our main results on the dynamical effect 
of accelerator modes on the diffusion exponent and the 
stable Levy distribution parameters. Linally, in Sec. 4 we 
discuss the main findings of this work. 

2. THE MODEL 

Lor the purposes of our study we use a rather well- 
known time-dependent system, the paradigm of classical 
chaos in Hamiltonian (nondissipative) systems, namely 
the kicked rotator [1], Its Hamiltonian function is given 
by 

H = + V 0 S T (t) cosO. (1) 

Here p is the (angular) momentum, I the moment of in¬ 
ertia, Vo is the strength of the periodic kicking, 9 is the 
(canonically conjugate, rotation) angle, and 5x(t) is the 
periodic Dirac delta function with period T. Taking into 
account that between the kicks the rotation is free, the 
Hamiltonian equations of motion can be easily integrated, 
and thus the dynamics can be reduced to the standard 
mapping, or so-called Chirikov-Taylor mapping, given by 

Pn+1 =Pn+V 0 sin0 n+ i, 

T (2) 

@n +1 — @n T —j^Pm 


Note that the quantities ( 9 n ,p n ) refer to their values just 
immediately after the n-th kick. Then, by introducing 
new dimensionless momentum P n = p n T //, we get 

P/i - l — Pn T K sin 9 n -\-i, 

dn+l = On + Pn, 

where the system is now governed by a single classical 
dimensionless kick parameter K = VqT/I, and the map¬ 
ping is area preserving. 

The generalized diffusion process of the standard map 
[Eq. (3)] is then defined as 

<(AP) 2 ) = D^K)n^ (4) 

where n is the number of iterations (kicks), and the expo¬ 
nent p is in the interval [0, 2), and all variables P, 9 and 
K are dimensionless. By D^{K) we define the general¬ 
ized classical diffusion constant and then according to the 
/;-value one may distinguish further the different types 
(rates) of diffusion. In more detail: in the case p = 1 
we have the normal diffusion , and DffK) is then the nor¬ 
mal diffusion constant, while in the case of anomalous 
diffusion we observe subdiffusion when 0 < p < 1 or 
superdiffusion if 1 < p < 2. In the case p = 2 we have 
the ballistic transport which is associated strictly - as we 
found - with the presence of accelerator modes. 

Whenever the diffusion is normal (p, = 1) the theoret¬ 
ical value of D\{K) is found to be [4] 

D (K) = f \K 2 [1 - *MK) (1 - MK))},if K > 4.5 
U |0.15(A' — A" cr ) 3 ,if K cr < K < 4.5 

(5) 

where K cr ,t ~ 0.9716 and ^(AT) is the Bessel function 
(the higher terms of order A' -2 are here neglected). 

The dependence of the diffusion constant for the 
growth of the variance of the momentum on K is not triv¬ 
ial and fails to be in a good agreement around the period 
1 accelerator mode intervals 

(27m) < K < (27m) 2 + 16, (6) 

where n is any positive integer (see [4, 5]). In more detail, 
in these intervals for the accelerator modes n = 1 there 
are two stable fixed points located at p = 0, 9 = n — 6 q 
and p = 0, 9 = n + 9 q, where 9q = arcsin(27r/AT). 
There are also two unstable fixed points atp = 0, 9 = 9q 
and p = 0, 9 = 2n — 9 0 . Furthermore, as the diffusion 
might even be anomalous, we have recalculated the effec¬ 
tive diffusion constant _D e g- = ((AP) 2 )/n numerically, 
which in general is not equal to the l) t , defined in Eq. (4) 
(more detail about this process can be found in [5]). 

3. RESULTS 

We have calculated the diffusion exponent p as a 
function of K. The final number of iterations is here 
n = 5000 and the total number of initial conditions is 
(a: 100000) spread uniformly on a 314 x 314 grid on the 
plane (0, P) = (0, 27t) (Fig. 1). The diffusion exponent p 
was calculated as described in Sec. II [Eq. (4)] in [5], i.e. 



Fig. 1 (Color online) The diffusion exponent p vs. K for 
an ensemble of psIOOOOO (314 x 314) initial condi¬ 
tions uniformly spread on the plane (0, P) = (0, 27 t). 
The intervals of stable accelerator modes of period 
1 are depicted on the black horizontal line p = 0.9 
while all the intervals of AT where the diffusion is nor¬ 
mal are found to have a diffusion exponent p k, 1). 
Note that the larger relatively peaks (for K > 27t 
- black full circles) correspond accelerator modes 
(mainly of period 1). Those peaks of smaller ampli¬ 
tude and for K < 27t (enlarged in the inset) as well as 
those for 2 tt < K < An marked with empty circle are 
related to accelerator modes of higher period. There 
are few K values marked with the symbol (x) asso¬ 
ciated with normal diffusion. All these cases where 
thoroughly studied in [5] where the interested reader 
may find more details. The power law describing the 
decay of the exponent p of the main peaks’ ampli¬ 
tude due to accelerator modes of period 1 (see text for 
more details) are shown by the blue dotted line. 


by the slopes of the lines of the variance of the momen¬ 
tum P as a function of iterations. The black horizontal 
line p = 0.9 indicates the intervals of stable accelerator 
modes of period 1 [Eq. (6)]. All intervals of AT asso¬ 
ciated with normal diffusion processes have p « 1. The 
large peaks, marked with full black circles, that appear for 
K > 27t are due to the anomalous diffusion raised by the 
accelerator modes of period 1. On the other hand, there 
are a few relatively smaller peaks for K < 2n (see the 
inset panel in Fig. 1), whose origin is accelerator modes 
of higher period, and also for 2n < K < 47t, both these 
sets are marked with empty circles. With the symbol (x) 
we mark few typical examples, close to those peaks, for 
which the diffusion is normal. All these cases where thor¬ 
oughly studied in [5] where the interested reader may find 
more details. 

In [5] it was found that the size of the large peaks 
for K > 2n (full black circles in Fig. 1 associated 
with regimes with accelerator modes of period 1) decays 
monotonically with a power law 

f(x) = cx d , (7) 

where c = 2.41645 and d = —0.195896 [blue dotted 




























line in Fig. 1, with asymptotic standard error ±0.04294 
(1.777%) and ±0.00537 (2.741%) respectively]. Hence, 
this shows that their effect decreases significantly for 
large A'-values, becoming probably entirely insignificant 
for K > 70. Moreover, the size of the successive acceler¬ 
ator modes of period 1 intervals decays with a power law 
defined simply and analytically by the Eq. ( 6 ). 

3.1. The Levy stable distribution 

Let us here give a brief description of the Levy stable 
distribution. For a random variable X with distribution 
function F{x), the characteristic function is defined by 

<j>{u) = Eexp(iuX) = J^° oo exp(iux)dF(x). Then, a 

random variable X is stable if and only if X = aZ ± b, 
with a > 0, b £ R and Z is a random variable with 
characteristic function 


E exp(kiZ) 


g—|u|“[l—i/3 tan ^(signn)] 
g—|u|[l+i/3tan j)(signu)] log |u| 


a ^ 1 
a = 1 ’ 
( 8 ) 


where 0 < a < 2 and — 1 < 0 < 1 (the symbol 
= indicates that both expressions have the same prob¬ 
ability law). We then adopt the parametrization k = 
0,«S'(a, 0, 7 , <5; 0) for which the random variable X given 
by 


Y ^ (7(^-/3tan(2^))±<5, a^l 
^7 Z + 6, a = 1 

has characteristic function 

S(a, 0, 7 , S; 0) = E exp(iuX) = 

( e iu<5-7 a |u|“(l-|-i/3(-l+|u7| 1- “)sign(u) tan(^)) 

\ -yM(-7r+2i/3 log(|tx 7 |))sign(M) 


(9) 


( 10 ) 

a ^ 1 
a = 1 


where Z = Z(a,f3) is defined as described in Eq. ( 8 ), 
a £ ( 0 , 2 ] is the index of stability or characteristic ex¬ 
ponent, [3 £ [— 1 , 1 ] the skewness parameter, 7 > 0 the 
scale parameter and S £ R location parameter. Two im¬ 
portant special cases are the Gaussian distribution with 
a = 2 and the Cauchy-Lorentz with a = 1 which are the 
only ones with an explicit closed formula. More details 
on the Levy stable distribution can be found in [ 6 ], 

Let us now see how one may utilize the Levy stable 
distribution in order to capture the effect of the anoma¬ 
lous diffusion due to the accelerator modes. In Fig. 2 
we depict the stable Levy distribution with parameters 
(a, 0 , 7 ,6) « (1.59,0.164,3.7,83.63) for A' = 6.5. In 
the inset color-map we show the exact grid of mixed ini¬ 
tial conditions (i.e. covering regular and chaotic regions) 
used (314 x 314 rH 00000). This particular cell contains: 
(i) trajectories that are transported ballistically by the un¬ 
stable accelerator mode around (0o,.P) = (1.31179,0) 
(x) and (ii) more evidently around the stable one at 
(0i, P) = (tt — 0o,O) (*) together with (iii) chaotic 
ones in their vicinity after n = 5000. The fit for the 
Levy distribution was performed by excluding the last 
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Fig. 2 (Color online) The distribution of the mo¬ 
menta P (black filled circles fitted by the Levy sta¬ 
ble distribution for K = 6.5. The best fit pa¬ 
rameters parameters are found to be (a, 0 , 7 , 6) ss 
(1.59,0.164,3.7,83.63) while the ensemble has 
~100000 (314 x 314) mixed initial conditions (shown 
in the inset), i.e. trajectories that are transported bal¬ 
listically by the unstable accelerator mode around 
( 00 , P) — (1-31179,0) (x) and more evidently 
around the stable one at (0i,P) = ( 7 T — 0o,O) (*), 
together with chaotic ones in their neighborhood af¬ 
ter n = 5000. In the inset color-map we show the 
diffusion exponent //, calculated by the process de¬ 
scribed in text, for a grid of cells on the subspace 
(0, P) £ [ 7 t/ 3 , 27t/ 3] x [ 0 , 7 r/ 3 ] of the phase space. 
In yellow (light gray in b/w) color we present ar¬ 
eas which correspond to ballistic motion due to the 
two accelerator modes pointed by the symbols (x) 
and (*) respectively. Those initial conditions being 
inside “pure” chaotic regime diffuse normally with 
p « 1. In the edges of the accelerator mode ar¬ 
eas where /i 0.8 one may find sticky subdiffusive 
transport (‘belt’’-like darker zone). The ballistic trans¬ 
port on the cylindrical phase space of the two above 
mentioned accelerator modes positioned initially at 
( 0 o,i, 0 ) and boosted by SP = 27r at every successive 
kick is shown in the small inset panel. 


peak in the positive large momenta due to the ballistic 
transport by the accelerator mode. Generally speaking, 
one should expect two peaks in the distribution of the 
diffusive variable P (in the positive and negative direc¬ 
tion) depending on the choice of the ensemble of ini¬ 
tial conditions. The color-plot inset panel shows the dif¬ 
fusion exponent /i for a grid of cells on the subspace 
(0,P) £ [ 7 t/ 3 , 2 - 7 t/ 3 ] x [ 0 , 7 r/ 3 ] of the phase space. The 
yellow (light gray in black and white) color areas are re¬ 
lated to ballistic motion due to the two accelerator modes 
pointed by the symbols (x) and (*) respectively. Those 
initial conditions lying inside the fully chaotic regime dif¬ 
fuse normally with p, « 1. Moreover, there is a “belt”- 
like darker zone in the edges of the accelerator mode 
area with p « 0.8 which indicates a “sticky” subdiffu¬ 
sive transport. In the small inset panel, one may see how 














the ballistic transport takes place on the cylindrical phase 
space of the two above mentioned accelerator modes po¬ 
sitioned initially at (0 o ,i, 0) and boosted by SP = 2n at 
every (only four here) successive kick. Furthermore, let 
us point out that the three yellow (or light gray in black 
and white) areas, surrounding the stable accelerator mode 
(*), are ballistically transported as well due to accelerator 
mode of higher period. 

Color maps like the one presented in the inset of Fig. 2 
can be used for the (numerical) detection of accelera¬ 
tor modes in the phase space. In Fig. 3 we show how 
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Fig. 3 (Color online) (a) the diffusion exponent /i for 
K = 3.1 and for 50 x 50 initial conditions on a 
500 x 500 cell grid of the entire phase space (0, P) £ 
[0, 27t] x [0, 27 t] calculated after n = 5000 iterations, 
(b) Two examples (marked with a square and circle in 
panel a) following different diffusion processes: a tra¬ 
jectory transported ballistically (with P < 0) by the 
effect of an accelerator mode of period 4 (circle) and 
one oscillating between islands of stability (square) of 
period 4 too. 

one can distinguish efficiently the regions with acceler¬ 
ator modes from the islands of stability with the use of 
diffusion exponent p color-plots. Let us first consider 
again a grid of 500 x 500 cells (on the entire phase space 
(i 0 , P) £ [0, 27t] x [0, 27 t]) with 50 x 50 initial conditions 
in each and evolve all of them together for n = 5000 
iterations. Then, for each ensemble of each cell sepa¬ 
rately, we calculate numerically the diffusion exponent p 
and append a color to this cell (color bar of Fig. 3(a)) that 
characterizes the different kind of diffusion of this small 
area. Depending on the location of each ensemble and 


the color assigned to the cell, one may expect to find: (i) 
normal diffusion (p = 1) inside chaotic regimes with¬ 
out the presence of accelerator modes, (ii) subdiffusion 
(0 < /i < 1) inside islands of stability, (iii) superdiffusion 
(1 < p < 2) inside chaotic regimes with the presence of 
accelerator modes in the phase space and (iv) ballistic 
transport (// ss 2) inside and in the very close vicinity of 
accelerator modes. 

Note that the stable regions around (7t/2,0), 
(37t/2, 0), (7t/ 2, 7r), (37 t/2, n r) and {it/2, 2n), (3n/2, 2t r) 
are correctly identified as islands of stability since their 
diffusion exponent p is smaller than 1. The remaining 
tiny yellow areas correspond to stable higher period ac¬ 
celerator modes with p « 2 (see [5] for more on the re¬ 
lation between the diffusion process and the stability of 
the different regimes of the phase space). Let us stress 
that following the above procedure, one manages to lo¬ 
cate the accelerator modes not only of period 1 but also 
of higher periods 2,3,4,... The latter ones are in general 
hard to be calculated analytically. Finally, in Fig. 3(b) 
we show two examples (marked with a square and a cir¬ 
cle in panel a) which present a rather different diffusion 
processes, i.e., for a trajectory oscillating between islands 
of stability (squares) and for one transported ballistically 
(with P < 0) by the effect of an accelerator mode (cir¬ 
cles). Let us notice here that the period of both is 4 (when 
projected on the (0, P) £ [0, 27t]), as it can be seen when 
we iterate them for a few steps. 

4. SUMMARY 

In this short paper we tried to give a small flavor of 
the accelerator modes’ impact in the diffusion process in 
the phase space of a system. These effects are not limited 
only to the standard map but also to other models that 
can have such dynamical features. One of the main re¬ 
sults is the link between the superdiffusive transport and 
the existence of the accelerator modes for certain control 
K parameter interval values. The interval size as well as 
the intensity (i.e. the deviation from normal diffusion ex¬ 
pressed by the deviation of the diffusion exponent p from 
p = 1) is found to decrease as a power law, as K in¬ 
creases (see Fig. 1). The effect of these modes seems to 
be strong even when their relative size in the phase space 
is small and even when one considers averages of initial 
conditions taken over the whole phase space. In the case 
where K > 27r the accelerator modes are in general of 
period one, while for K < 2tt they may be of higher pe¬ 
riod, and their influence on p decays with K much faster 
than for period one accelerator modes. 

Looking into these resonances of the phase space even 
further, we tried to use an adequate probability distribu¬ 
tion function of the momenta P that would allow us to 
distinguish between normal and anomalous diffusion for 
different ensembles of initial conditions. It turns out that 
the stable Levy distribution and its a parameter can serve 
for this goal quite well, despite the fact that this distri¬ 
bution is mainly introduced and broadly used to capture 
random (Levy) walks in the literature. Here, we showed 






only one exemplary case for K = 6.5 (see Fig. 2) where 
the tt-Levy parameter was found to be w 1.59 while the 
/j diffusion exponent lies on a range of values depending 
on the location of the sub-cell of initial conditions (inset 
of Fig. 2). 

Finally, by considering a fine grid in the phase space 
and then taking ensembles of many initial conditions for 
the averages inside the small cells, we managed to pro¬ 
duce diffusion “chart-map” encoding their different (nu¬ 
merically estimated) // value as shown for example in 
Fig. 3. With this kind of plots one may detect and sep¬ 
arate accelerator modes (p ~ 2) of different periods from 
even small islands of stability (0 < p < 1), or chaotic 
regions (// « 1), or “sticky” intermediate regimes. 

We should mention that another important method to 
study the structure of the phase portrait rests upon cal¬ 
culating the Lyapunov exponents in order to distinguish 
the regular from the chaotic regions in the phase space, 
and also allows us to quantify the degree of chaos by the 
value of the Lyapunov exponents [7] at the sufficiently 
densely and uniformly chosen initial conditions. How¬ 
ever, the calculation of the Lyapunov exponents can be 
quite laborious task, while the so-called GALI method, 
a generalization of the original SALI method, introduced 
and developed by Skokos et al [8-11], and recently re¬ 
viewed by Skokos and Manos in reference [12], is an ex¬ 
tremely efficient method, which distinguishes quickly the 
regular from chaotic orbits, although it does not quantify 
the degree of chaos to such an extent as the Lyapunov ex¬ 
ponents. This method has been used in our main work 

[5] and it was found that there is excellent agreement 
between the //-landscape and the GALI-land sc ape in the 
phase space of the standard map. We believe that our ap¬ 
proach and the method can be applied in other smooth 
Hamiltonian dynamical systems and mappings. 
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